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POTENTIAL USE OF HIGH-FREQUENCY INDUCTION HEATING FOR HIGH- 
TEMPERATURE LIQUID-METAL HEAT-TRANSFER TESTING 
by John V. Miller 
Lewis Research Center 

SUMMARY 

Experimental determination of the local heat -transfer coefficients of liquid metals at 
temperatures in excess of 1000° K is very difficult. Not only must the heaters used in 
the test operate at a temperature where structural strength is low, but also, because of 
the corrosiveness of liquid metals at these temperatures, the test section must be made 
of a refractory metal and, therefore, be opeictieu m a comroneu tjnvixunmeiu 10 prevent 
oxidation. 

To investigate potential methods of obtaining high-temperature liquid-metal heat- 
transfer data, several possible heating devices have been examined and a cursory sum- 
mary of the problems of each are discussed. One method which appears promising is 
the use of high-frequency induction heating. The advantage of this type of heater is that 
the heating coil is physically separated from the test section and can be maintained at a 
much lower temperature than the test section. This allows the ceil to be constructed of 
relatively low -temperature materials. 

An analysis of a thick-walled cylinder, internally cooled with flowing liquid metal 
was performed to evaluate the characteristics of an inductively heated test section. It 
was assumed that a suitable heating coil could be constructed for this application. The 
results of this study indicate that, if the ratio of cylinder radius to current depth is 
greater than 6, internal heat generation in the center of the cylinder is negligible. With 
a tungsten cylinder, for example, commercially available high-frequency equipment can 
be used to obtain the necessary cylinder-radius to current-depth ratio without exceeding 
allowable thermal stress levels which have been established for tungsten. 

Neglecting the temperature variation of electrical resistivity in the analysis does not 
appreciably influence the calculated temperature distribution of a tungsten cylinder since 
nearly all the induced power is in a region where the temperature and resistivity varia- 
tions are small. Based on an error analysis of an inductively heated test section, ex- 
perimental determination of liquid-metal heat-transfer coefficients at high heat -flux 
levels can be made to an accuracy of ±10 percent for reasonable tolerances on thermo- 
couple placement and temperature measurements. An additional error is introduced if 
the thermal conductivity of the material used in the test section is not well defined. 


INTRODUCTION 


Liquid metals are excellent heat-transfer fluids because of their relatively high 
thermal conductivity, high volumetric heat capacity and low vapor pressure at elevated 
temperatures. These characteristics make liquid metals ideal candidates for many high- 
temperature, high heat -flux applications. In particular, the use of liquid metals as 
coolants for nuclear reactors is being considered extensively for systems requiring 
coolant temperatures of 1000° to 2000° K. In this temperature range, lithium, sodium, 
potassium, rubidium, cesium, and a mixture of sodium and potassium (NaK) appear to 
have the necessary qualifications for nuclear reactor coolants (refs. 1 to 3). 

In liquid-metal cooled systems, a large portion of the heat transferred is a result of 
the high molecular conductivity of the coolant. Therefore, heat-transfer correlations 
for ordinary fluids do not apply since with most other fluids the molecular conductivity is 
negligible compared with the eddy conductivity. 

Much work has already been done in the field of liquid-metal heat transfer. However, 
the majority of the experimental work has been performed at relatively low temperatures 
(<1000° K) and low heat fluxes (<200 W/cm^). Extrapolation to higher temperatures and 
higher heat fluxes should be possible except that there is wide scatter in the available 
data. Figure 1, for example, shows some of the many heat -transfer relations which 
have evolved from theoretical (refs. 4 to 6) and experimental (refs. 7 and 8) studies of 
various liquid metals. Although the form of these and many similar liquid-metal heat- 
transfer equations vary considerably, they all fall roughly into the region between the 
Martinelli and Kulateladze equations shown on figure 1. 

Generally, the maximum variation in the predicted heat-transfer coefficient is about 
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a factor of 2 or varies by ±35 percent about an average line which could be constructed 
for the various relations (e. g. , at a Peclet number of 1000, the predicted heat -transfer 
coefficient for lithium flowing in a 5-ram tube at 1000° K is about 12±4 W/cm 2 -°K). For 
many applications, this ±35 percent accuracy is sufficient. At a heat flux of 100 watts 
per square centimeter and a nominal heat -transfer coefficient of 5 watts per square 
centimeter per °K, for example, the difference between the fluid temperature and pre- 
dicted surface temperature is 20°±7° K. The relatively small error (±7° K) could be 
considered negligible for cases in which such temperature differences were not critical. 
For advanced nuclear reactors, however, heat fluxes in excess of 400 watts per square 
centimeter and possibly as high as 800 to 1000 watts per square centimeter are desirable 
to increase the specific power output and/or to minimize the size and weight of a given 
reactor system, particularly for space-power applications where weight is at a premium. 
At these high heat fluxes, a ±35 percent error in predicting a nominal heat-transfer co- 
efficient of 5 watts per square centimeter per °K could result in an error of more than 
50° K which, under certain conditions, could result in a reactor fuel element exceeding 
some limiting heat -transfer or material strength criteria. For example, if a lithium 
cooled reactor were operating at a coolant temperature of 1750° K and at a pressure of 
3 atmospheres (3. 03x10^ N/'m 2 ), a 50° K error in predicting the surface temperature 
could produce a film boiling burnout condition (ref. 9) since the surface temperature 
could be higher than the saturation temperature of lithium at that pressure (i. e. , 

>1790° K). To prevent such adverse situations, it appears that additional liquid-metal 
heat -transfer experiments will be required to refine the present data if very high specific 
power reactors are to be developed. 

The reasons for the wide scatter in presently available liquid-metal heat -transfer 
data are varied. 

(1) It is difficult to obtain accurate heat -transfer data at high temperatures and heat 
fluxes because of instrumentation problems. 

(2) For liquid-metal heat-transfer tests, the system must be clean and free from 
oxides and other contaminants. Much of the scatter in the data has been attributed to the 
fact that this condition has not always been achieved (ref. 10). 

(3) Another important factor, which is of special concern at the high temperatures 
of interest for nuclear reactors, is that the high corrosivity of liquid-metals dictates the 
use of refractory metals (i. e. , niobium, molybdenum, tantalum, or tungsten) and some 
of their alloys as the material of construction. Since the refractory metals are suscep- 
tible to oxidation, they must be used in a controlled environment (inert gas or vacuum) 
further increasing an already difficult instrumentation problem. 

The instrumentation problem is not merely related to the ability to accurately mea- 
sure temperatures under high flux and temperature conditions, but it also involves selec- 
tion of the best locations for measuring temperatures in the test and the accuracy of 
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thermocouple placement. This factor has several aspects which are peculiar to the var- 
ious methods used to obtain heat-transfer data in liquid-metal systems. The purpose of 
this report is to briefly discuss some of these methods and to explore the potential use 
of high-frequency induction heating for liquid-metal heat -transfer testing. 


METHODS OF OBTAINING LIQUID-METAL HEAT-TRANSFER DATA 

The task of obtaining reliable experimental heat-transfer data for liquid metals at 
high temperatures is more difficult as the temperature level increases. At the temper- 
atures of interest for space-power nuclear reactors (1000° to 2000° K), even the selec- 
tion of the best method of testing is difficult. Because the liquid metals are to be tested 
at a high temperature (e. g. , 2000° K), the test section and heater must be at a signifi- 
cantly higher temperature in order to obtain the required 400- to 1000-watts per square 
centimeter heat flux. 

Although there are several methods of testing which can possibly be applied to 
liquid-metal heat transfer, the high-temperature problem limits the applicability of 
many of these methods and reduces the effectiveness of others. Some of these testing 
methods and their related problems are discussed in the following sections. 


Resistance Heating 

Probably one of the most widely used methods of obtaining heat-transfer data is the 
use of direct resistance heating of the test section using direct -current or low-frequency 
alternating- current voltage. In such a system, the coolant is circulated through a rela- 
tively thin-walled tube, a section of which serves as both heater and test section. Ther- 
mocouples are attached to the outer surface. The heat-transfer coefficient between tube 
and coolant is determined from the temperature measurements and a heat balance of the 
system . 

Obviously, since the thermocouples are on the outside of the tube, the temperature 
of the inside (heat transfer) surface must be estimated. This could lead to a significant 
error in the predicted heat-transfer rate, particularly for large heat fluxes and tempera- 
ture gradients. However, for applications in which the fluid being tested is not a good 
electrical conductor, this method has been used successfully to determine the heat- 
transfer characteristics of gases (ref. 11) and of liquids such as water (ref. 12). But, 
when the fluid being studied is a very good electrical conductor, such as is the case of 
liquid metals, the current being applied is split between the flowing fluid and the test 
section. A portion of the total heat load is, therefore, generated directly in the coolant. 
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It would be necessary to correct the temperatures and surface heat flux accordingly to 
accurately establish the heat-transfer coefficient. Because the determination of the 
electrical current split between test section and fluid requires knowledge of the electrical 
resistivity of both materials and because the resistivities are generally temperature de- 
pendent, an accurate correction is, in practice, very difficult. 


Insulated Resistance Heating 

One method of overcoming the current shunting problem would be to separate the 
liquid-metal coolant from the electrically heated test section by use of an electrical insu- 
lator. Several such systems have been designed and tested. One of them uses a 0. 3- 
centimeter -diameter tantalum heater, insulated from a 0. 95-centimeter-diameter 
niobium -zirconium sheath by 0. 18 centimeter of swaged alumina (ref. 13). This heater 
has operated 500 hours at a surface flux of 20 watts per square centimeter with a center- 
line temperature of approximately 1700° K. Even with this relatively thick insulator 
some electrical leakage occurred; the effective resistivity of the alumina was reduced to 
approximately 25 percent of the handbook value. 

A second type of resistance heater for liquid-metal systems (ref. 14) consists of a 
graphite core, approximately 5 centimeters long and 0. 3 centimeter in diameter, which 
is contained in a boron nitride insulator and sheathed in stainless steel. Several of these 
heaters have operated at fluxes in excess of 470 watts per square centimeter for over 
400 hours. The use of stainless steel, of course, limits the allowable coolant tempera- 
ture to less than 1000° K. 

Although further advances in size, heat flux, and allowable operating temperature of 
these insulated resistance heaters will undoubtedly be achieved, there is one basic factor 
which limits their application to very high temperatures and heat fluxes. This is the 
relatively large thermal gradient that will occur across the insulator at high flux levels. 
For example, at a heat flux of only 500 watts per square centimeter, the temperature 
drop across a 0. 2-centimeter alumina insulator is about 1250° K. If the liquid-metal 
coolant were operating, for example, at 2000° K, it is obvious that the heater temperature 
(>3250° K) would be well over the melting point of even the refractory metals. There- 
fore, the use of an insulated electrical heater for high-temperature, high heat -flux appli- 
cations is limited by the allowable temperature of the heater and, to an extent, by the 
possible breakdown of the electrical insulating properties of the insulator at high temper- 
atures. 
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Thick-Wall Resistance Heating 


Another way to minimize the current shunting problem without using an insulator is 
to use a thick-walled rather than a thin-walled tube for the heater and test section. The 
electrical resistance of any element is proportional to the resistivity and length, and in- 
versely proportional to the cross-sectional area (appendix B). For the case in which the 
resistivity of the liquid-metal coolant and the thick-walled test section have about the 
same value, the fraction of heat generated directly in the liquid metal is approximately 
equal to the ratio of the areas of the coolant passage and the cylinder (i. e. , (q^/q^.) = 
(a/R)^ ). Figure 2 illustrates how the fraction of heat generated directly in the liquid 
metal decreases as the thickness of the tube is increased relative to the coolant hole size. 
When the diameter of the test section is 10 or more times the size of the coolant hole 
(i.e., R/a > 10), over 99 percent of the heat is generated in the tube walls, and the heat 
generated in the liquid-metal coolant could be considered negligible. 

In this illustration (fig. 2), it was assumed that the resistivity of the tube was the 
same as that of the liquid metal. In practice, such a condition would probably not occur, 



Figure 2. -Effect of cyl in der-to -coolant -hole size on heat generated in coolant 
of resistively heated, liquid-metal -cooled test section assuming that electrical 
resistivity of liquid metal and the test section are approximately equal. 
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since, even within the thick-walled cylinder, the resistivity would vary with temperature. 
This illustration, however, does point out that, by increasing the wall thickness suffi- 
ciently, the shunting problem can be neglected. 

Because of the thick-walled tube used in this method, it would be impractical to at- 
tach thermocouples to the outer surface and make accurate corrections for the tempera- 
ture on the inside surface. This is apparent from figure 3 where the temperature distri- 
bution in a resistively heated cylinder shows a significant difference in the surface tem- 
perature for typical conditions which might be encountered in a test. A difference of 
nearly 70° K exists at the outer surface of the cylinder between the case where all heat 
generated in the test section is transferred to the fluid (e = 1.0) and the case where 
50 percent of the heat is lost by conduction, convection, and radiation from the outer 
surface. 

A more logical method of instrumenting the thick -walled tube would be to drill a hole 
radially from the outer surface to as near the inner surface as possible. The thermo- 
couple could be then positioned near the inner wall in a location where the thermal effi- 
ciency e of the test section does not affect the temperature profile. A heat balance in 
the liquid-metal coolant could establish both the fluid temperature and heat flux at the 
same axial location as the thermocouple. This information would suffice to establish the 
local heat -transfer coefficient. A correction would, of course, be necessary to account 
for the difference between the actual inside surface temperature and the temperature at 



Figure 3. - Typical temperature distributions in resistively 
heated thick-walled cylinder cooled by liquid metal. 

Ratio of cylinder to coolant-hole size, 12; radius of 
coolant channel, 0.32 centimeter; bulk fluid tempera- 
ture, 1670° K ; heat-transfer coefficient, 6 watts per 
square centimeter per °K; heat flux, 475 watts per 
square centimeter. 
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the thermocouple location. A discussion of the accuracy with which the heat-transfer 
coefficient can be measured is given in the section ERROR ANALYSIS. 


DOUBLE TUBE HEAT EXCHANGER 

Average heat-transfer rates can be obtained using a double-tube heat exchanger. In 
this method an inner tube serves as a carrier for one fluid, while a second fluid passes 
through the annulus formed by the two tubes. Heat transfer from the second fluid to the 
inner tube can be by single-phase convection or by condensation of a vapor on the outer 
wall of the inner tube. The secondary fluid may or may not be the same as the primary 
fluid in the first case; it may actually be the identical fluid which is pumped in series 
through both the outer and inner tube. In the case where the fluids are different, the 
secondary fluid may even be a gas, although it would require a substantial flow of gas to 
provide the necessary heat -transfer coefficient and heat capacity to supply the required 
heat load to the liquid metal. 

The chief problem with this system is that it is impractical to attach thermocouples 
on the inner tube. Temperatures measured on the outer surface of the outer tube would 
need considerable correction to obtain a local heat -transfer coefficient for the inner tube 


Thermal-Radiation Heaters 

The thermal-radiation heater is one in which the heat transfer to the test section is 
by thermal radiation from a much hotter heat source. One such heater has been designed 
(ref. 13) to produce about 20 watts per square centimeter at 1000° K. 

The main disadvantage of this type of heater to high-temperature, high heat -flux ap- 
plications is the very high operating temperature of the heater. For example, at a cool- 
ant temperature of about 1500° K, the required heater temperature for a thermal- 
radiation device would be over 3500° K to deliver 1000 watts per square centimeter using 
the ideal black-body equation (ref. 15): 

9 = F°(tJ - T<) (1) 

(Symbols are defined in appendix A. ) 
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Electron- Bombardment Heating 


An electron-bombardment heater consisting of an evacuated tube in which a stream 
of electrons were directed against the inner wall has been used (ref. 13) with liquid- 
metal systems. This heater had already operated at 110 watts per square centimeter but 
only has an expected capability of 300 watts per square centimeter. 

One difficulty with electron-bombardment heating is the complexity necessary to en- 
sure a uniform stream of electrons over the entire surface of the test section. For ex- 
ample, to obtain uniform heating on the outside of a round tube it would be necessary to 
use multiple electron beams to effectively cover the entire surface. Further work on 
electron-bombardment heaters may eliminate this restriction and improve the perfor- 
mance of this type of heater. 


High-Frequency Induction Heating 

Another potential method of performing high-temperature heat -transfer tests on 
liquid metals is the use of high-frequency (>1000 Hz) induction heating. Induction heating 
is widely used industrially on many configurations for such applications as brazing, an- 
nealing, tempering, and many other processes in which localized high-temperature heat- 
ing effects are required. The rapidly alternating current of such a system results in the 
so called skin-effect which produces a highly concentrated heating source near the sur- 
face of the metal piece. The power input could either be directly coupled across the test 
section, as is presently done with direct-current or low-frequency alternating current 
systems, or the more conventional method of inducing a current in the test section by 
passing a current through a surrounding coil could be used. Because most of the power 
is generated near the surface, the heat transfer near the center of the test section is es- 
sentially a conduction process. 

Kalachev, et al. (ref. 16), have investigated the use of induction heating for flowing 
liquid metals which are to be used as heat-transfer fluids. They have concluded that, 
for high-temperature application, such a heater has the advantages of smoothly adjust- 
able output, stable heat generation, and a power output which is independent of the tem- 
perature of the heated fluid. They have estimated the heat flux capability for this appli- 
cation at about 1000 watts per square centimeter and have recommended induction heating 
for use in all liquid-metal systems. 

However, Kalachev, et al. , were mainly interested in using the induction heater as a 
means of heating the liquid metal. To extend this work to the experimental determination 
of the local heat -transfer coefficients of liquid metals, it is first necessary to determine 
the basic thermal and electrical characteristics of an inductively heated test section used 
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Pump 



Figure 4. * Diagram of liquid-metal heat -transfer loop using high-frequency induction heating. 


for such an application. The resulting temperature distribution and related information 
obtained from an analysis of an inductively heated test section can be used to define the 
optimum location for placement of thermocouples to enhance the reliability of the experi- 
ment and to reduce the error in determining the heat -transfer coefficient. 

Figure 4 schematically illustrates the major components required for an inductively 
heated system. The coolant loop consists of pump, flow control valve, and a heat ex- 
changer. By regulating the flow rate in the primary coolant loop and in the secondary 
side of the heat exchanger, the liquid metal could be maintained at some predetermined 
temperature and flow level for each test. 

The instrumented test section is a thick-walled cylinder (fig. 5) surrounded by the 
high-frequency induction coil. To avoid disturbing the induced electrical field, the 
thermocouples would probably have to be inserted axially and located near the center of 
the cylinders in the region where the induced current flow is negligible. This could pos- 
sibly result in some error due to inaccuracies in locating the thermocouple (see ERROR 
ANALYSIS). 

The advantage of this type of heater is that the induction coil may be physically sepa- 
rated from the test section by insulation or some form of radiation shield to reduce heat 
losses from the hot test section and to allow the coil to operate at a much lower tempera- 
ture than the test section. Induction coils are generally constructed of hollow copper 
tubing, which is water cooled to limit the operating temperature of the copper. 
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Figure 5. - Test section. 


However, because of the temperature difference between test section and the sur- 
rounding coil, a portion of the heat generated in the test section is lost from the outer 
surface of the cylinder. The ratio of the heat actually transferred to the liquid metal to 
that generated in the test section will be termed the thermal efficiency e of the test sec- 
tion. In addition to the heat loss directly from the test section, other losses occur in the 
system. One source of such losses is the coil itself. Although induction heating is used 
quite extensively in many industrial applications, in many cases the actual construction 
of the coil is an art rather than a true science. To obtain uniform heating and to prevent 
excessive end losses, the coil must be delicately tuned for each application. However, 
for a cylindrical test section the construction of a suitable coil does not appear to be a 
major problem, although the operating frequency must be properly selected to avoid a 
mismatch between the impedance of the power supply and the test section. 

Kalachev, et al. (ref. 16), have measured the overall efficiencies of inductively 
heated liquid-metal systems at 55 to 75 percent. With such efficiencies, the power re- 
quirement for the high-frequency generator would be nearly twice that actually needed 
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for the test. For example, to transfer 1000 watts per square centimeter to a 25- 
centimeter -long test section with a 0. 5-centimeter coolant-hole diameter, only 39 kilo- 
watts of power would actually be transferred to the liquid metal, but the required power 
output of the high-frequency generator would be about 70 kilowatts (i. e. , using Kalachev's 
55 percent overall efficiency). However, this is well within the power levels of commer- 
cially available equipment as shown on table I. 


TABLE I. - COMMERCIALLY AVAILABLE HIGH-FREQUENCY 
POWER SUPPLIES 


Power supply 

Frequency range, 
kHz 

(a) 

Power output range, 
kW 

(a) 

Refer- 

ence 

Electronic tube generators 

180 to 450 

1 to 100 

17 

Electronic tube generators 

2500 to 25 000 

2. 5 to 50 

17 

Motor -generators 

1 to 10 

7.5 to 500 

18 


a Not all combinations of frequency and power are available in the ranges 
shown. 


Therefore, from the standpoint of being able to obtain the necessary power for a 
liquid-metal heat -transfer test, it appears that presently available commercial equipment 
could be utilized. It remains necessary to show how well the heat-transfer coefficient of 
liquid metals can be determined from an inductively heated test section. 


ANALYSIS OF INDUCTIVELY HEATED TEST SECTION 
Induced Current 


The fundamental equation for the distribution of the magnetic -field intensity H in an 
inductively heated cylinder is (ref. 19) 


£H + ISH-k 2 JH = 0 


dr‘ 


r dr 


( 2 ) 


where k is defined by the relation 


k 2 _ 2vfjJ. 

P( T ) 


(3) 
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The induced current is related to the magnetic -field intensity by 


I 


= J_ dH 
4 7r dr 


( 4 ) 


Since nearly all nonferrous metals have about the same permeability (ref. 20) 

(p n = 47 tX 10”9 Wb/A-cm) a thick-walled refractory metal cylinder with a liquid-metal 
coolant will behave essentially the same as a solid cylinder with respect to the induced 
current distribution. This is particularly true when nearly all the induced current is 
near the outer surface of the cylinder, in which case the presence or absence of the cool- 
ant in the center of the cylinder is unimportant to the current distribution. 

The solution to equations (2) and (4) is complicated by the variation of the electrical 
resistivity p of most metals with temperature. Table n shows the resistivity of some 

TABLE n. - RESISTIVITY OF SOME TYPICAL 
REFRACTORY METALS 3 


Temperature, 

°K 

Tantalum 

Molybdenum 

Tungsten 

Resistivity, microhm -cm 

293 


5.7 

5.51 

300 

13.85 

— 


1000 

44. 1 

— 

25.3 

1023 


21.0 


1500 | 


— 

41.4 

2000 


— 

59.4 

3000 


— 

98.9 

3510 


— 

118.0 


a Data from ref. 21. 


typical refractory metals (data from ref. 21). In the temperature range from 1500° to 
3000° K, for example, the resistivity of tungsten can be represented by an equation of the 
form 


P( T ) 


2.333 


T_V 


1000/ 


+ 27.83/-l-\- 5.587 

\iooo/ 


10 ohm -cm 


(5) 


Obviously, substitution of this expression into equation (2) would result in a very complex 
differential equation. And, since the functional form of the temperature dependency of 
resistivity is different for each metal, a completely general solution to the inductively 
heated cylinder would not be possible. To obtain the current distribution in an inductively 
heated cylinder, therefore, the assumption is made that the resistivity of the cylinder is 
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independent of temperature. The validity of this assumption can be determined after the 
complete solution for the inductively heated cylinder, including the resulting temperature 
profile, has been obtained. 

Reference 19 shows that, under the assumption of uniform resistivity, the ratio of 
the induced current at any radial location r to that induced at the surface of the cylinder 
is given by 


*(r) _ ber'(kr) + j bei'(kr) ^ 

I ber’(kR) + j bei’(kR) 

where ber and bei are defined, respectively, as the real and imaginary parts of the 
Bessel function of the first kind of zero order; that is, 


J Q (jkr\/j"^ = ber(kr) + j bei(kr) (7) 

Values of ber and bei and the derivatives ber' and bei' are tabulated in various 
forms in references 19, 22, and 23. The induced current distribution may be numeri- 
cally evaluated using these tables and considering the absolute value of equation (6a); 
that is, 



ber ,2 (kr) + bei ,2 (kr) 
ber ,2 (kR) + bei ,2 (kR) 


1/2 


(6b) 


Heat Conduction 

If the length-to-diameter ratio of the cylinder is large (L » d), the temperature 
gradient in the axial direction will be small compared with the radial gradient. Under 
these conditions, axial heat conduction away from the ends of the cylinder can be ne- 
glected, and the radial, steady- state heat conduction in the cylinder with internal heat 
generation (ref. 24) is given by 

d 2 T , 1 dT . Q (r, T) (8 

dr 2 r dr k (T) 

As in the case of the induced current, the temperature dependency of the properties 
of the material complicate the solution of the differential equation. Therefore, the as- 
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sumption is made that the resistivity and thermal conductivity are independent of temper- 
ature. The validity of this assumption will also be discussed after the temperature pro- 
file has been determined. 

Appendix B shows that, for a uniform resistivity, the local volumetric heat- 
generation rate is related to the induced current by 

%) . [Vf 

Qo Jo_ 


so that 


Q(r) Q 0 


ber t2 (kr) + bei ,2 (kr) 
ber ,2 (kR) + bei ,2 (kR)_ 


( 9 ) 


It is obvious that the direct use of equation (9) in the heat conduction equation (eq. (8)) 
would also result in a rather complex differential equation. 


Approximate Solution 

Reference 19 notes that for large values of the quantity kR, an approximate form of 
equation (6a) is given by the exponential expression 


(r) _ -k(R-r)/Vi 

1 — V 


(10a) 


When the value of R - r is equal to V2/k in this expression (eq. (10a)), the local value 
of the induced current is reduced by a factor of e or to 36. 8 percent the value at the 
outer surface of the cylinder; that is, when (R - r) = -y/I/k, 


l (r) _ 1 _ 


= ± = 0.368 
e 


For this reason, the value of (yfe/k) has been termed (ref. 19) the equivalent cur- 
rent depth 6 of an inductively heated cylinder and is defined as the depth, relative to the 
outer surface of the cylinder, at which the induced current is reduced by a factor e. 
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TABLE IE. - VARIATION IN INDUCED CURRENT AND VOLUMETRIC HEAT 


GENERATION RATE OF INDUCTIVELY HEATED CYLINDER WITH 
INCREMENTAL VALUES OF EQUIVALENT CURRENT DEPTH 


Incremental value of equiva- 
lent current depth a 

Relative value of induced 
current' 3 

Relative value of heat 
generation rate 13 

1 

0. 3679 

1.354X10" 1 

2 

.1353 

1.83 xlO' 2 

3 

.0498 

2.5 xlO" 3 

4 

.0183 

3.0 xlO' 4 

5 

.0067 

<10' 4 

6 

.0025 

<10' 5 

7 

.0009 

<10' 6 

8 

.0003 

<10" 7 


a Measured from outer surface of cylinder. 
^Relative to the value at outer surface of cylinder. 


(For each successive incremental 6 distance, the current is reduced by a factor of e 11 
as shown on table m. ) 

Equation (10a) can be rewritten in terms of this defined variable 6. 


V) = e -(R-r)/6 


where 



(10b) 


(11a) 


which, for the case of nonferrous materials (j u. = 4rrxl0 


-9 


Wb/A-cm), 


can be expressed as 


6 = 5033.7 



(lib) 


A plot of the variation in induced current with cylinder radius is shown on figure 6 for 
several values of the dimensionless ratio R/6. It can be seen that, for the larger values 
of R/6, the induced current is reduced to a negligible value in a relatively short distance. 

Because the volumetric heat -gene ration rate is related to the square of the induced 
current (eq. (9)), the falloff in heat generation is even more rapid (table HI). Figure 7 
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Dimensionless radius, r/R 

Figure 6. - Variation in induced current with radius and 
cylinder radius to current depth ratio using approximate 
solution (eq. (10b)). 



0 .2 .4 .6 .8 1.0 

Dimensionless radius, r/R 


Figure 7. - Comparison between exact and approximate 
solution of volumetric heat-generation distribution in 
cylindrical body located in high-frequency induction 
coil. 



0 .2 .4 .6 .8 1.0 

Dimensionless radius, r/R 


Figure 8. - Relation between exact and approximate volu- 
metric heat -generation distribution for cylindrical 
distribution for cylindrical body located in high- 
frequency induction coil. 
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illustrates this point and also compares the volumetric heat -gene rat ion distribution ob- 
tained using the exact solution (eq. (9)) and that obtained using the approximate form of 
the current distribution (eq. (12)): 

5 (£) = e - 2 ( R - r >/8 ( 12 ) 

% 

For the larger values of R/5, the absolute difference between the exact and approxi- 
mate solutions is small. Equally important, the exact solution for a given value of R/6 
lies between the approximate solution for that and another value of R/6. This can be 
seen more readily in figure 8, where the exact solution for kR = 10 lies between the 
approximate solutions of kR = 8 and kR = 10. This is an important factor since, as 
will be shown, the resulting temperature distributions for either value (kR = 8 or 
kR = 10) is essentially identical in certain regions of the cylinder. Therefore, since 
these two approximate solutions bracket the exact solution for kR = 10, the temperature 
distribution for the exact solution would also be identical in these regions. 


TEMPERATURE DISTRIBUTION 

Using the approximate solution for the volumetric heat generation (eq. (12)), the 
solution for the radial temperature distribution in an inductively heated cylinder is de- 
rived in appendix B. The result, in terms of the dimensionless temperature difference 

K[T/ x - tJ 

AT/ v = -LM (1! 

j a <?a 


is given by 
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where 


x = 2- 
5 

y = 2 - 
5 

and e is the thermal efficiency of test section (i.e. , the ratio of heat actually entering 
the coolant to total heat generated in test section). 

Figure 9 shows the radial temperature distribution which results from applying the 
dimensionless expression (eq. (14)) to a set of conditions which might be typical of an in- 
ductively heated test section. Figure 9(a) shows the variation in temperature for several 
values of the dimensionless ratio R/5. Also shown in figure 9(a) is the temperature 
distribution which would result if the heat transfer from the outer surface of the cylinder 
to the coolant were only by conduction (i. e. , no internal heat generation). Note that, for 
larger values of R/6, the temperature distribution in the inductively heated cylinder ap- 



(a) For several values of ratio of cylinder radius to cur- (b) For test section thermal efficiencies of 0.50 and 1.00. 

rent depth. Thermal efficiency, 0.50. Ratio of cylinder radius to current depth, 6. 

Figure 9. - Radial temperature distribution in inductively heated cylinder. Ratio of cylinder radius to coolant hole size, 12. 
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proaches that of the conduction distribution. Conversely, as the value of R/6 decreases, 
the temperature distribution becomes more like that of the thick-walled resistively heated 
cylinder (fig. 3). 

These temperature profiles (fig. 9(a)) were determined using a test section efficiency 
of 50 percent; that is, half of the heat generated in the test section actually enters the 
coolant and half is lost from the outer surface of the cylinder. Although the power re- 
quirement of the test section is inversely proportional to this efficiency, figure 9(b) shows 
that, for large values of R/5, the effect of test section efficiency on the temperature dis- 
tribution of an inductively heated cylinder is small. In fact, in the region r/R < 0. 6, 
there is essentially no difference among the temperature distributions of (1) an inductively 
heated cylinder with a test section thermal efficiency of 0. 5, (2) an inductively heated cyl- 
inder with a test section thermal efficiency of 1.0, or (3) a cylinder in which heat transfer 
is by conduction only. The reason for this similarity is shown on figure 10 where lines of 
constant power -input for an inductively heated cylinder are shown as a function of cylinder 
radius and the ratio R/6. For values of R/6 > 6, over 99 percent of the total heat is 
generated in the region 0. 6 < r/R <1.0. Therefore, for a given heat load to the coolant, 
the temperature distribution in the region r/R <0.6 is essentially that of a conduction 
process and is independent of conditions near the outer surface of the cylinder. This be- 
havior is an important factor in the application of induction heating to heat -transfer test- 
ing. 



0 .2 .4 .6 .8 1.0 

Dimensionless radius, r/R 


Figure 10. - Lines of constant power output for inductively heated 
cylinder. Percentage of heat generated between a given radius 
and outer surface of cylinder. 
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APPLICATION TO HEAT TRANSFER TESTING 


In applying induction heating to tests in which local heat -transfer coefficients are to 
be determined, it is necessary to establish 

(1) The heat flux 

(2) The surface temperature of the coolant passage 

(3) The coolant temperature 

Ideally, for purposes of correlating experimental data, it would be desirable to as- 
sume that the heat generation within the cylinder is negligible and that the heat transfer 
from the outer surface of the cylinder to the coolant channel was by conduction only. As 
already demonstrated (figs. 7 to 10), this is a fair assumption provided that the ratio of 
cylinder radius to current depth R/6 is properly chosen. For values of R/6 greater 
than 6, the temperature distribution in the region a/R < r/R <0.6 is essentially the 
same as that which would result from heat conduction through a thick-walled cylinder 
with no internal heat generation. 

The heat flux can be determined by measuring the temperature difference between 
two thermocouples located in this region (a/R < r/R < 0. 6) and applying the heat conduc- 
tion equation (ref. 24) for an internally cooled cylinder, 


K(T 2 - Tj) 


\n( T J\ 

W 


(15) 


where the subscripts 1 and 2 refer to thermocouple locations (fig. 5). 

When equation (15) is rewritten in terms of the coolant -channel surface, the temper- 
ature at the inner surface of the thick-walled cylinder can be determined using either of 
the two thermocouples as a reference 


T 


a 



<?a 



K 


(16) 


The coolant temperature at the axial location of interest can be measured directly 
or, can be established by measuring the inlet- and exit-coolant temperature and assuming 
a linear temperature rise across the test section. (Note that this method can be used in 
conjunction with the flow rate to check the heat flux. ) 
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Once these three parameters cp„, T„, and T, have been established, the heat- 

a a D 

transfer coefficient can be calculated from the relation 

w 

h = (17) 

< T a - T b> 

The only restriction to using the above method of obtaining heat -transfer data for an 
inductively heated cylinder is the requirement that the R/5 ratio be sufficiently large to 
ensure that the heat -gene rat ion rate in the region of the thermocouples is negligible. 

There are two ways in which the ratio R/6 can be altered to obtain the required value 
(R/6 > 6): either the diameter of the cylinder can be increased, or the equivalent cur- 
rent depth can be decreased. Obviously, the advantage of using a large diameter cylin- 
der for the test section is that it provides additional cross-sectional area to facilitate the 
placement of thermocouples (fig. 5). However, the larger the diameter of the cylinder, 
the greater will be the temperature gradient through the test section for a given power 
input to the coolant. For example, for a 5-, 7. 5-, and 10-centimeter tungsten cylinder, 
the temperature gradients associated with a 500-watt -per-square-centimeter heat flux at 
the surface of a 0. 75 -centimeter cooling hole are 435°, 580°, and 685° K, respectively. 
Although it is not clear that such gradients could not be tolerated, the thermal stresses 
imposed, particularly for cyclic operation, could result in test section failure. 

Fiero (ref. 25) has established allowable temperature gradients for tungsten based 
on cyclic strain data. It appears that, for 25 to 50 cycles, temperature gradients of ap- 
proximately 500° K could be tolerated without test section failure. Therefore, depending 
on heat flux and coolant-hole diameter, tungsten cylinders of about 6 to 8 centimeters in 
diameter could be utilized for an inductively heated test section. To maintain a radius- 
to-current depth ratio R/5 of 6 or more for cylinders of this size, it would be neces- 
sary to keep the equivalent current depth below 0. 5 to 0. 67 centimeter. Figure 11 shows 
a plot of required frequency necessary to obtain various current depths. Over the tem- 
perature range of interest (1000° to 2000° K), the minimum frequency necessary to main- 
tain the required current depth (5 < 0. 67 cm) is 1500 hertz. At frequencies above this 
value, the equivalent current depth is decreased, and the resulting R/5 ratio is further 
improved (i. e. , increased to values greater than 6). Since this frequency range (i. e. , 
>1500 Hz) is available on commercial equipment (table I), inductively heated tungsten 
cylinders, 6 to 8 centimeters in diameter, appear feasible for a liquid-metal heat- 
transfer test section. Table IV lists the characteristics of a typical tungsten test section 
for an arbitrary set of conditions. 
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TABLE IV. - TYPICAL CHARACTERISTICS 21 OF AN INDUCTIVELY 
HEATED TUNGSTEN TEST SECTION FOR LIQUID -METAL 
HEAT-TRANSFER TESTING 


(a) General characteristics 


Coolant hole diameter, cm 

0.64 

Test section length, cm 

25 

Heat -transfer area, cm^ 

50 

Outside diameter of test section, cm 

7.68 

Test section-to -coolant hole ratio 

12 

Minimum cylinder radius-to-current depth ratio 

6 

Maximum allowable current depth, cm 

0.64 

Required operating frequency, Hz 

b 3500 

Average liquid -metal temperature, °K 

1670 

Heat flux to coolant, W/cm 2 

475 

Heat -transfer coefficient, W/(cm 2 )(°K) 

6 

Surface temperature of coolant channel, °K 

1750 

Effective thermal conductivity of test section, W/(cm)(°K) 

b l. 03 


(b) Power requirements as function of thermal efficiency 


Power requirements 

Thermal efficiency of 
test section, percent 


100 

50 

Maximum test section temperature, °K 

2100 

2090 

Power required by test section, kW 

23.75 

47.50 

Minimum power supply requirement, c kW 

31 

62 


a For arbitrary set of conditions. 

^Assuming an average test section temperature of 2000° K. 
c Assuming a 30-percent loss in power supply and heating coil. 


ERROR ANALYSIS 

The local heat-transfer coefficient of an inductively heated, liquid-metal cooled cyl- 
inder can be determined using the methods outlined in the previous section. However, 
the accuracy with which the heat-transfer coefficient can be established depends on sev- 
eral factors. 


Temperature Variation of Properties 

The first item which could influence the accuracy of determining the heat -transfer 
coefficient, and the most important to this analysis, is the extent to which the assumption 
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Figure 12. -Typical temperature distributions in induc- 
tively heated tungsten, thick -walled cylinder convec- 
tively cooled by liquid metal. Ratio of cylinder radius 
to current depth, 6; ratio of cylinder radius to coolant- 
hole size, 12; radius of coolant channel, 0.32 centi- 
meters; thermal efficiency, 0.50; bulk fluid tempera- 
ture, 1670° K ; heat-transfer coefficient, 6 watts per 
square centimeter per °K. 

of constant resistivity and constant thermal conductivity can be applied. 

Figure 12 shows typical temperature distributions in an inductively heated tungsten 
cylinder for conditions which might be encountered in a liquid-metal heat-transfer test. 
These temperature distributions were determined using the equation which was derived 
with the assumption of constant properties (eq. (14)). The variation in electrical resis- 
tivity of tungsten with temperature has already been given (eq. (5)). A similar equation 
(ref. 26), expressing the thermal conductivity of tungsten in the temperature range of 
interest, is given by 


K/,^ = 133. 8 - 15. 6 -I- W/(m)(°K) (18) 

1000 


Although the slope of the temperature distribution in the inductively heated cylinder 
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Resistivity, p, (ohmMcm) Thermal conductivity, 
K, (W/(m)(°K)> 
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(a) For several ratios of cylinder radius to current depth. 
Ratio of cylinder radius to coolant-hole size, 12; heat 
flux, 475 watts per square centimeter; thermal effi- 
ciency, 0.50. 
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(b) For test section thermal efficiencies of 0. 50 and 1. 0. 
Ratio of cylinder radius to current depth, 6; ratio of 
cylinder radius to coolant-hole size, 12; heat flux, 

475 watts per square centimeter. 
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(c) For several heat-flux levels. Ratio of cylinder radius 
to current depth, 6; ratio of cylinder radius to coolant- 
hole size, 12; thermal efficiency, 0. 50. 

Figure 13. - Variation in resistivity and thermal conductivity with test section radius. Radius of coolant channel, 0. 32 centi- 
meter; bulk fluid temperature, 1670° K; heat -transfer coefficient, 6 watt per square centimeter per °K. 
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appears to be quite steep, particularly for high values of heat flux (fig. 12), the resulting 
variation in the properties is not. Figure 13 shows the variation is thermal conductivity 
and electrical resistivity across the cylinder for various values of R/5 , thermal effi- 
ciency e and heat flux. The extreme variation in thermal conductivity for the worst 
case shown (i.e. , for the highest heat flux) is less than 10 percent; for lower values of 
heat flux the variation is even less. 

In the case of the electrical resistivity, the gross variation across the cylinder is on 
the order of 40 percent for the high heat flux case (fig. 13(c)). However, in the region 
where most of the heat is generated and where the variation in resistivity is important, 
the variation about a mean value of resistivity is only ±2 percent for the 790-watt -per- 
square -centimeter heat -flux case and ±1 percent for the 475-watt -per-square -centimeter 
case under the conditions shown. (Note that even if the heat -transfer coefficient were 
changed by a factor of 2 from the 6 W/(cm^)(°K) used in this example, the absolute tem- 
perature level would only change by a small amount because the temperature difference 

between surface and coolant is small compared with the temperature drop through the 

o 

tungsten. At 475 W/cm , for example, the temperature level would only change by 
40° K for a factor of 2 increase in the heat -transfer coefficient. ) Obviously, the higher 
the R/5 ratio and the lower the heat flux, the less will be the variation in resistivity in 
the power producing region. 

It can be concluded that neglecting the temperature variation in properties should 
not have an appreciable effect on the analysis of a tungsten cylinder for the range of 
variables covered in this section. However, each material and each set of conditions 
must be checked to ensure that the property variation in the region of interest is small. 


Negligible Internal Heat Generation 

Another factor which could influence the accuracy of a heat-transfer test, using the 
method outlined in the section Application to Heat Transfer Testing, is the assumption 
that the heat generation in the region r/R <0.6 is negligible. The validity of this as- 
sumption has qualitatively been established for large values of R/6. A quantitative ex- 
pression can also be obtained. 

The heat conduction equation for the case of no internal heat generation in a thick- 
walled cylinder is given by equation (16) which can be rearranged and expressed in terms 
of the dimensionless temperature ratio (eq. (13)) as 

AT =ln- (19) 

c a 
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Comparison of equations (14) and (19) shows that the dimensionless temperature differ- 
ence of the inductively heated cylinder (eq. (14)) is equal to that of the value of conduc- 
tion case (eq. (19)) plus a correction term. Defining this correction term as a yields 


a = (ATj - AT C ) 



The ratio of the induction to conduction temperature distributions, which is a measure of 



0 .2 .4 .6 .8 1.0 0 .2 .4 .6 .8 1.0 


Dimensionless radius, r/R 

(a) Test section thermal efficiency, 0. 50. <b> Ratio of cylinder radius to current depth, 6. 

Figure 14. - Effect of current depth on difference in radial temperature distribution between inductively heated and conduc- 
tively heated cylinders. Ratio of cylinder radius to coolant-hole size, 1?. 
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the importance of the internal heat generation, can be then written as 


y = i + — (21) 

ml 

a 

\ 

Figure 14 shows plots of the dimensionless temperature ratio y for several values 
of R/5 and e . For values of R/6 greater than 6, the difference between the tempera- 
ture distribution for the inductively heated cylinder and the conduction case are negligible 
in the region r/R < 0. 6, and the internal heat generation can therefore be neglected in 
this region. 


Experimental Errors 

In addition to the aforementioned possible sources of error which could result from 
assumptions made in the analysis, there are also errors which could result from experi- 
mental inaccuracies. These errors include those resulting from measuring tempera- 
tures and temperature differences and those resulting from insufficient knowledge of 
thermocouple locations. Experimental errors are common to all heat-transfer tests of 
this nature and are not peculiar to the induction heating method. However, the fact that 
the thermocouples must be inserted axially (fig. 5) in the inductively heated test section 
could result in a somewhat larger error than if they were inserted radially. 

The magnitude of the experimental errors can be determined by combining equations 
(15) to (17) into a dimensionless heat transfer coefficient 



In addition to the desired heat-transfer coefficient, the left side of this equation con- 
tains the thermal conductivity of the test section and the radius of the coolant hole. Al- 
though these two parameters directly influence the accuracy of the experiment, they are 
essentially a part of the test section rather than the experimental measurements. The 
error resulting from these parameters will largely depend on how well the properties of 
the test section are known. 

Grouped on the right side of the dimensionless relation (eq. (22)) are those factors 
which directly affect the experimental data. These include the location of thermocouples 
and the ability to measure temperatures and temperature differences. Figure 15(a) 
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-.6 -.4 -.2 0 .2 .4 .6 -.6 -.4 -.2 0 .2 .4 .6 

Error in thermocouple 2 placement, mm 

(a) Location of thermocouple 1, 1. 15 centimeters. (b) Location of thermocouple 1, 0. 4 centimeter. 

Figure 15. - Effect of error in thermocouple location on ability to predict heat -transfer coefficient from measured temperature difference. Radius 
of coolant channel, 0. 32 centimeter; ratio of cylinder radius to current depth, 6; bulk fluid temperature, 1670° K; heat flux, 475 watts per square 
centimeter; thermal conductivity, 1.03 watts per square centimeter per °K. 

shows how an error in thermocouple placement causes an error in the predicted value of 
(ha/K). These calculations were performed for an initial thermocouple placement (1) of 
r i = 1. 15 centimeters and for various separation distances (rg - r ^). The more distance 
between the thermocouple locations, the less sensitive the data is to an error in thermo- 
couple placement. Figure 15(b) shows that, by placing the thermocouples nearer the 
coolant hole (i. e. , r^ = 0. 4 cm instead of 1. 15 cm), the error is reduced considerably; 
for an error in thermocouple placement of -0. 6 millimeter, for example, the error in 
predicting ha/K for a nominal separation distance of 1.0 centimeter is reduced from 
13. 5 to 5 percent for the conditions shown. 

The effect of error in measuring temperature differences is shown on figure 16. Al- 
though the slopes are reversed, the magnitude of the error involved is essentially the 
same for either an error in (Tg - Tj) or (Tg - T^). Since there are many possible com- 
binations of errors in temperature measurements and thermocouple locations, it is im- 
possible to estimate the exact error which would result in a typical experiment. Table V, 
however, shows the resulting error for several heat-flux levels and for all possible com- 
binations of these errors assuming that the temperature differences could be measured to 
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Error in measured temperature difference, °K 


Figure 16. - Effect of error in measuring temperature on the ability 
to predict heat-transfer coefficient. Radius of coolant channel, 

0. 32 centimeter; location of thermocouple 1, 0.4 centimeter; dis- 
tance between thermocouples 1 and 2; 1.2 centimeters; bulk fluid 
temperature, 1670° K; heat flux, 475 watts per square centimeter; 
thermal conductivity, 1.03 watts per square centimeter per °K. 


±3° K and that the location of the thermocouple r 2 “ r j = !• 2 centimeters could be deter- 
mined to ±0. 6 millimeter. The range of errors for a heat flux of 790 watts per square 
centimeter is about ±10 percent, while that for 475 watts per square centimeter is 2 to 
4 percent higher. 

Based on these values, it would appear that an overall error of ±20 percent could be 
achieved if the thermal conductivity of the test section were reasonably established (e. g. , 
to ±10 percent). 


COMPARISON WITH THICK-WALLED RESISTANCE HEATING 

Basically, the experimental errors involved in the thick-walled resistively heated 
test section are the same as those for the inductively heated case (figs. 15 to 17). How- 
ever, an additional complication arises because of the temperature variation of the elec- 
trical resistivity. In the inductively heated case, the temperature variation in the region 
of maximum heat generation was small enough that the resistivity variation was negligi- 
ble (fig. 13(c)). However, in the resistively heated test section, the heat is generated 
throughout the cylinder and the entire variation in the resistivity must be considered. 
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Figure 17. - Comparison of resistivity variation in tung- 
sten cylinder heated by induction or resistance heating 
and cooled by liquid metal. Radius of coolant channel, 

0. 32 centimeter; ratio of cylinder radius to coolant-hole 
size, heat flux, 475 watts per square centimeter; 
bulk fluid temperature, 1670° K ; heat transfer coeffi- 
cient, 6 watts per square centimeter per °K. 

Figure 17 illustrates this point; the resistivity variation for an inductively and a resis- 
tiveiy heated thick-walled cylinder are compared for a set of assumed operating condi- 
tions. 

For the resistance heater, the variation in the electrical resistivity across the cyl- 
inder is between 18 and 23 percent depending on the thermal efficiency of the test section, 
hi the inductively heated case, the variation in resistivity is less than three percent for 
the region in which 99 percent of the heat is generated. It is obvious that the result of 
neglecting resistivity variation will be much more significant for the resistively heated 
test section. 


CONCLUDING REMARKS 

Based on the analysis of an inductively heated, thick-walled cylinder, induction heat- 
ing appears to offer a potential method of obtaining high temperatures (1000° to 2000° K) 

O 

and high heat fluxes (400 to 1000 W/cm ) for liquid-metal heat -transfer experiments. 

The analysis shows that, if the ratio of cylinder radius to equivalent current depth is 
greater than 6 and if the thermocouples are located in the region (a/R < r/R < 0. 6), in- 
ternal heat generation may be neglected. The resulting temperature distribution in this 
region is essentially the same as if the heat transfer through the cylinder were by con- 
duction only. 

Although no attempt was made to actually design a heating coil for the test section 
discussed in this report, it was assumed that a practical coil could be constructed which 
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would operate with a reasonable efficiency (e. g. , 50 percent). Under these conditions, 
commercially available high-frequency equipment will adequately supply the power re- 
quirements for high heat-flux testing and provide the frequency necessary to achieve the 
desired cylinder radius to current depth ratios using a tungsten cylinder approximately 
6 to 8 centimeters in diameter. 

Assuming that the variation in electrical resistivity in a tungsten test section may be 
neglected does not affect the results of the analysis because the variation in temperature 
and resistivity is small in the region where the heat generation occurs. However, this 
assumption must be checked for a given set of conditions if materials other than tungsten 
are used or if the heat flux is very high. 

Errors in the determination of the heat -transfer coefficient which are introduced by 
inaccurate thermocouple placement and temperature measurements are on the order of 
10 to 15 percent over the heat flux range of interest (400 to 1000 W/cm ). An additional 
error can also be introduced if the thermal conductivity of the test section is not well 
established. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 2, 1967, 

120-27-06-17-22. 
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APPENDIX A 


SYMBOLS 


a 

radius of coolant channel, cm 

P 

electrical power, W 

ber,bei 

defined parameters of Jq 

Pe 

Peclet number, GDeC p /K 

C 1’ C 2 

constants of integration (ap- 
pendix B) 

Q 

volumetric heat -generation 
rate, W/cm 3 

C P 

specific heat at constant pres- 
sure, W-sec/(g)(°K) 

% 

heat loss from outer surface of 
cylinder, W 

D 

mathematical operator, d/dz 

q 

heat generated, W 

De 

equivalent hydraulic diameter, 

R 

outer radius of cylinder, cm 


cm 

R* 

dimensionless distance, r/a 

d 

diameter of coolant channel, 

(ft 

electrical resistance, ohms 


cm 

r 

arbitrary radius, cm 

F 

geometric view factor for 
radiant heat transfer 

Ar 

tolerance on thermocouple 
location, mm 

f 

frequency, Hz 

T 

temperature, °K 

G 

9 

mass flow rate, gm/(sec)(cm ) 

AT 

dimensionless temperature 
difference, k[t^ - T a ]/ a< P a 

dimensionless maximum tem- 

H 

magnetic field intensity, A/cm 

AT* 

h 

I 

heat transfer coefficient, 
W/(cm 2 )(°K) 

electrical current, A 


perature difference, 
K < T max - Ta’/^a 



X 

dimensionless parameter, 

J 0 

Bessel function, first kind- 


2 (a/5) 


order zero 

y 

dimensionless parameter, 

i 

imaginary unit, V-l 


2(R/5) 

K 

thermal conductivity, 
W/(cm)(°K) 

a 

dimensionless temperature 
difference, AT^ - AT c 

k 

variable defined by eq. (3) 

y 

dimensionless temperature 

L 

length of test section, cm 


ratio, AT i /AT c 

Nu 

Nusselt number, hDe/K 

5 

equivalent current depth, cm 

n 

integer 1, 2, 3, 4, . . . 

e 

thermal efficiency of test sec- 
tion 
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M 

permeability, Wb/A-cm 

max 

maximum value 

P 

resistivity, (ohm)(cm) 

n 

nonferrous materials 

cr 

Stefan-Boltzmann constant, 

0 

outer surface of test section 


5. 7X10" 12 W/(cm 2 )(°K 4 ) 

o-r 

region from center of cylinder 

<P 

heat flux, W/cm 2 


to r 

Subscripts: 

r 

radius (variable) 

a 

value at coolant surface 

r-R 

region between r and R 

b 

bulk fluid condition 

T 

temperature 

c 

conduction only 

t 

total value 

h 

heater 

1,2 

thermocouple locations 1 and 2 

i 

induction heating 
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APPENDIX B 


DERIVATION OF EQUATIONS FOR INDUCTION HEATING OF CYLINDRICAL 
TEST SECTION FOR HIGH-TEMPERATURE HEAT-TRANSFER 
TESTING OF LIQUID METALS 

Equation for Temperature Distribution in Cylindrical Body 
in Inductive Electrical Field 


If it is assumed that the electrical resistivity of the material is constant, the approx- 
imate solution for the induced current distribution in a cylindrical body (ref. 19) is given 
by 


^W. e -«-r)/S (bu 

: o 

for large values of R/5, where 1^ is the current at any radial location and I Q the 
value at the surface of the cylinder. (The result of neglecting the resistivity variation is 
discussed in the main text. ) 

The power dissipated by the flow of current in a resistive element is 

P = I 2 <ft (B2) 

where the electrical resistance (ft is proportional to the resistivity of the material and 
the length of the resistor and inversely proportional to the cross-sectional area of the 
resistor 


(ft = p (length/area) (B3) 

For the cylindrical shell (fig. 5), the geometry of the resistive element is 

area = dr L 
length = 2irr 
volume = 27rrL dr 
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so that 


<r =iH£ 

L dr 


and 


p _ 27rrp t 2 
P (D - ^ >(r) 

The volumetric heat -generation rate of the cylindrical shell is given by 


Q 


(r) 


P (r) _ P (r) 


volume 2?rrL dr 


Combining equations (B5) and (B6) 


Q (r) = 


2?rrpI (r) 

L dr 27rrL dr 


Q (r) 


(L dr)‘ 


l (r) 


(B4) 


(B5) 


(B6) 


(B7) 


Again neglecting temperature variations in the resistivity of the material, the volumetric 
heat-generation rate may be related to the induced current 


Q (r) _ 

Q o 



(B8) 


or 


Q( r ) = Q o e" 2(R ' r)/6 (B9) 

Neglecting axial conduction and temperature variations in thermal conductivity, the basic 
heat conduction equation for a cylindrical body with internal heat generation is given by 
reference 24. 



1 dT 
r dr 


Q (r) 

K 


(BIO) 
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which, for the case of induction heating, becomes 



1 dT 
r dr 


-Q e -2(R-r)/5 

n) 

K 


letting 


A = e -2R/6 

K 


and 



equation (Bll) may be expressed as 


d 2 T 


dr 


A Al = Ae mr 
r dr 


or 


2 d 2 T 
r - + 

dr 2 


_ dT . 2 mr 
r — = Ar e 

dr 


This has the general form of the Cauchy linear differential equation which can be 
by letting 


r = 


e 


z 


or 


z = In r 


Equation (B12) may then be rewritten as (ref. 27) 

D(D - 1)T + DT = Ae 2z e meZ 

where D is the mathematical operator d/dz. Expanding equation (B13) yields 


(Bll) 


(B12a) 


(B12b) 

solved 


(B13) 
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or 


(D 2 - D + D)T = Ae 2 z e meZ 


D 2 T = Ae 2z e me 


= Ae 2z e meZ 

dz 2 


The complementary solution of this equation is 


T = CjZ + Cg 


To obtain a particular solution, let 


u = 


dT 

dz 


so that 


and 


Letting 


and 


= Ae 2z e me 
dz 


"- A / 


e 2z e me dz 


s = e 


ds = e dz 


equation (B16) can be written as 


u = A J' se ms ds 


(B14) 


(B15) 


(B16) 
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for which 


Therefore, 


and 


Again letting 


ms 

u = A (ms - 1) 

m 2 

m _ 

r = /A^y e me2 (me z - l)dz 

s = e z 
ds = e z dz 


equation (bT 7) becomes 



so that 



By substituting for s, the particular solution is found to be 



(B17) 


(B18) 
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Combining the complimentary (eq. (B15)) and particular solution (eq. (B18)), the com- 
plete solution is given by 


T (z) “ c i z + c 2 + 


\m 


e me -Z 


/ j n -n! 
n=l 


replacing Z by In r, equation (B19) becomes 


T ( r) = c i lnr + C 2 + | 

The previously substituted values of 




mr 


- tor -V 
A n ‘ n! 


n=l 


Q o -2R/6 


and 


A = - — e 
K 


m 2 
m = — 

6 


may be resubstituted at this time so that 


T W - Cl in, + C, (|\ 2 


In r - e 


2r/6 + 



2 - 


n 

6 

n -n! 


n=l 


The first derivative of this expression is also useful in the solution of the problem 

m-l' 


dT = ^1 + 5° e -2R/5/5 
dr r k \2 


— - e 2r / 6 + 
2r 



2r 

6 


n! 


n=l 


(B19) 


(B20) 


(B21) 


(B22) 
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The boundary conditions of interest in the problem of a liquid metal heat -transfer test 
are 


and 


<p a = K — at r = a 
a dr 


T a = T (r) at r = a 


When these boundary conditions are substituted into equations (B21) and (B22), the tem- 
perature distribution in the test section is 


T (r)- T a = C l 




e- 2R / 6 «!ln0 + e 2a / 6 


V 



(B23) 


where 


C 


1 




e -2R/6 



Determination of Heat Generation Rate at Cylinder Surface 

In the derivation of the temperature distribution in the cylinder, it was assumed that 
the value of Q q , the heat generation rate at the outer surface of the cylinder, was known. 
However, the value of Q q is dependent on several factors and must be related to these 
variables to obtain a complete solution to the problem. 

The total heat generated in the cylinder by induction heating can be found by inte- 
grating the volumetric heating rate 


r R 

q t = y o 2^rLQ (r) dr 


(B24) 


43 


which on substituting the value of and integrating becomes 


q t = 2 ttLQ o 0 2 ^e _2R/6 + 2^-1^ (B25) 

In the same manner, the integrated heat input between any two arbitrary radial locations 
may be determined. The results for the regions of interest are 


q 


o-r 


= 27tLQ 0 


/6\ 2 e -2R/5 

~ e 2r/8 / g £ . A , 

\2 / 

\ 8 ) J 


so that 


q 


o-r 


* 


e 2 r/6 

e 2R/6 



+ 1 


+ 1 


(B26) 


(B27) 


and 


1.0-5“^ (B28) 

q, q, 


Since some heat may be lost from the outer surface of the cylinder by radiation, conduc- 
tion, and convection, the total heat generated in the cylinder is the sum of the heat lost 
from the outer surface and that transferred to the liquid-metal coolant 

q^ = 27raL(p a + q^ (B29) 

Defining the thermal efficiency of the test section as the ratio of heat received by the 
liquid-metal coolant to the total heat generated in the cylinder, 


2naL<p 

a 


(B30) 


(This is different than the overall efficiency of the system where losses in other compo- 
nents must also be considered. ) Using this relation, equation (B29) can be written as 
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(B31a) 


^t = 


27iaL<p a 

e 


or 


2ttL e 


(B31b) 


Substituting in equation (B25), the expression for Q q is found in terms of the heat flux to 
the coolant 


a* e 2R / 6 



(B32) 


Dimensionless Temperature Distribution 

When equation (B32) is substituted into equation (B23), the complete solution for the 
temperature distribution is obtained 



(B33) 


where 
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Defining the dimensionless temperature difference 

K[T/ * - T 1 

AT = L (r) M (B34) 

a< ?a 

and rearranging equation (B33) yields an expression for the temperature distribution 
which is independent of the heat flux to the coolant and the thermal conductivity of the 
material. 



(B35) 


where the dimensionless parameters x and y are given by 


and 


x 




46 


Determination of Maximum Temperature in Cylinder 

The maximum temperature T max in the cylinder will occur when 

d(AT) _ Q 



By differentiating equation (B35), a transcendental expression for the location R* of the 
maximum temperature difference AT* is obtained 


XR*e XR 


00 



n=l 



(B36) 


where R* = (r/a) at the location of maximum temperature. This equation may be solved 
graphically or by trial and error and the value of R* substituted into equation (B33) to 
obtain the expression for the maximum temperature difference 


r 


AT* = ln(R*) + lJ 


e x [l+X ln(R*)J - e xR * + 



e y (y - 1) + 1 



(B37) 


v. 


J 
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